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Abstract: We study SO{d+ 1) invariant solutions of the classical vacuum Einstein 
equations in p + d + 3 dimensions. In the limit d —)■ oo with p held hxed we con¬ 
struct a class of solutions labelled by the shape of a membrane (the event horizon), 
together with a ‘velocity’ held that lives on this membrane. We demonstrate that 
our metrics can be corrected to nonsingular solutions at hrst sub-leading order in 
i if and only if the membrane shape and ‘velocity’ held obey equations of motion 
which we determine. These equations dehne a well posed initial value problem for 
the membrane shape and this ‘velocity’ and so completely determine the dynamics 
of the black hole. They may be viewed as governing the non-linear dynamics of the 
light quasi normal modes of Emparan, Suzuki and Tanabe. 
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1 Introduction 


The rich classical dynamics of uncharged black holes is governed by the vacuum 
Einstein equation 

R^u = 0 . ( 1 . 1 ) 

This innocuous looking equation captures very complicated processes; for example 
the collision of two black holes and their subsequent merger accompanied by grav¬ 
itational radiation. We owe much of our current understanding of these processes 
to numerics. However the numerics involved are very challenging; their complexity 
makes it impractical to densely hll the space of initial conditions with numerical solu¬ 
tions. It would thus be useful to have analytic techniques to analyse these solutions. 
However the phenomena described are of such complexity that it seems unlikely that 
exact analytic solutions will ever be obtained. Moreover perturbation theory appears 
to be ruled out by the fact that Einstein’s equations lack a parameter. 

In this paper we follow the lead of Emparan, Suzuki and Tanabe (EST) to intro¬ 
duce a parameter into Einstein’s equations. We do this by studying these equations 
in D dimensions. The equations of black hole dynamics simplify in this limit and 
admit a systematic expansion in Our general strategy is similar to that of t’Hooft 
[1] who introduced a parameter into the study of Yang Mills theory by replacing the 
SU{3) gauge group by SU{N), and that of Witten who introduced a parameter into 
the study of atomic and molecular physics by analysing the quantum mechanics of 
the ^ potential in large dimensions [ 2 ]. 

EST and collaborators have recently pointed out [3-10] that several features of 
black hole dynamics simplify in the limit of a large number of spacetime dimen¬ 
sions D. A Schwarzschild black hole with Schwarzschild radius tq in D spacetime 
dimensions is described by the metric 


ds'^ 






( 1 . 2 ) 


If r is held hxed at a value greater than tq as H —)■ cx) then ^ ^ 0 and the 

metric (1.2) reduces to flat space. On the other hand if we set r = ro(l -|- 
and hold R hxed as D is taken to inhnity then ^ —)■ e~^. It follows that 

the gravitational tail of a black hole in D spacetime dimensions extends only over 
a distance away from its event horizon [3], a distance we will refer to as the 
‘thickness of the membrane’. 

EST have computed spectrum of quasinormal modes of the black hole (1.2) 
in an expansion in A [ 5 ^ 3 ]. They hnd that almost all of the inhnite number of 
quasinormal modes at every angular momentum have frequencies of order the inverse 
of the membrane thickness, i.e. of order R. However at each angular momentum 
a small number of modes (two scalars and one vector at each value of the angular 
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momentum) are much lighter than the generic mode; the frequencies of the light 
modes are of order Moreover the light quasinormal modes are supported entirely 
inside the membrane region. ^ 

The separation in scale between the heavy and light quasinormal modes suggests 
that the full nonlinear dynamics of a black hole over time scales much larger than 
^ is governed by an effective non-linear theory of the light quasinormal modes. We 
expect the degrees of freedom of this theory to reside within the membrane region, 
and so effectively on a codimension one surface in spacetime at length scales large 
compared to 

In this paper we will verify the expectation of the previous paragraph by explicit 
analysis of Einstein’s equations. We identify the auxiliary dynamical system that 
lives on the membrane and determine its dynamical equations of motion. Our work 
has several similarities to the Fluid-Gravity correspondence ([11] see [12, 13] for 
reviews), however there are also signihcant differences. The construction of this 
paper applies to gravity even in flat space. The effective ‘hydrodynamical’ theory 
obtained in this paper lives on a fluctuating surface rather than the hxed held theory 
background of huid gravity. Finally the constructions presented in this paper are [14] 
justihed not by an expansion in gradients in units of the horizon radius, but rather 
by an expansion in T. Our membranes also has several conceptual similarities to the 
blackfolds (see e.g. [15],[16], [14], [17]), but differ from the later in decoupling from 
radiation. 

A solution of Einstein’s equations in D dimensions is given by a metric tensor 
with D{D + l)/2 components; each of these components is a function of D spacetime 
variables. It is clear that the generic solution of this sort has no analogue at any hxed 
hnite value of D. For example a metric in four dimensions has many fewer legs and 
and many fewer directions in which to wiggle than a metric in a million dimensions. 
In order to take a sensible large D limit we write D = p + d + 3 and divide up the 
D spacetime dimensions into two groups. The hrst group contains p + 2 dimensions 
including time and the second contains the remaining d + 1 dimensions. We restrict 
our attention only to solutions that preserve the SO{d + 1) rotational symmetry in 
d + 1 dimensions. We then take the limit d ^ oo with p held hxed. ^ 

In the rest of this paper we use the following notation. Let Wa represent the hrst 
set of p -|- 2 coordinates. Let zm denote the last set of d -|- 1 coordinates. The metric 
of hat D dimensional spacetime is 

= dwadw°' + dzMdz^ = dwadw°‘ + dS'^ + 

^In contrast the heavy quasinonrmal modes are non-trivial outside the membrane region where 
they reduce to purely outgoing modes. 

^It follows that the solutions we study have analogues in all dimensions greater than or equal to 
p -t 3. 
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Let = {S, w°') so that the index jj, runs over p + 3 variables. The most general 
metric that preserves the SO{d + 1) isometry takes the form 

ds'juii = g^vdx^dx'" + (1.4) 

where and (j) respectively are a metric and a scalar held in the p + 3 dimensional 
space spanned by in particular they are functions only of x^ and not of the 
coordinates of the d sphere. 

Below we derive the equations of motion for g^y and 0, and explain that these 
equations admit an interesting large d limit. We then explicitly construct a large 
family of solutions of these equations. The solutions we construct are in one to one 
correspondence with the solutions of an auxiliary problem of membrane dynamics in 
hat space that we now describe. 

Our auxiliary dynamical system lives on a codimension one membrane that re¬ 
sides in the hat space (1.3). This membrane is free to huctuate, subject to the 
requirement that it preserves SO{d+ 1) invariance. The position of the membrane 
can always be characterized by the zeroes of a function B{x^) and SO{d -|- 1) in¬ 
variance forces B{x^) to be a function only of and S (and not of the angular 
coordinates on the d sphere). We could, for example, choose the function B to take 
the form B = S — g{w°') so that the membrane surface is given simply by the equa¬ 
tion S = g{w°‘). It follows that the data in the shape of the membrane surface is 
contained in a single function of p -|- 2 variables. 

For many purposes we can think of the membrane as propagating in the auxiliary 
hat space 

ds^ = dwadw°‘ + dS"^ (1.5) 

Let denote the outward pointing unit normalized normal one form held of the 
membrane, is assumed always to be spacelike, is of course completely deter¬ 
mined by the shape of the membrane surface; in terms of the function B we have 

h* y/dB.dB ■ 

Our membrane is also equipped with a ‘velocity’ held on its world volume 
is a p -|- 3 dimensional one form held that is further constrained by the following 

^However the S coordinate is inequivalent to the other coordinates and so our effective 
dynamical system enjoys invariance only under the p + 2 dimensional Lorentz invariance. 

^More precisely our membrane has a null oneform field O living on its world volume. It turns out 
that the overall scale of O can be absorbed into a subleading shift of the membrane surface and so 
is irrelevant at leading order. Upto a scaling this oneform can be cast in the form O = dS — Uadw’^; 
the fact that 0.0 = 0 then implies that u.u = —1. This normalization leads us to refer to the 
collective coordinate as a ‘velocity’. The precise physical significance of Uy, will be clear once we 
understand its contribution to the membrane stress tensor. We leave this to future work. 
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unusual conditions 


u^ = -l 

u.n = ns - (1.6) 

ns 

Us = 0 

(here ns and us are simply the S components of the one forms n and u). The hrst of 
( 1 . 6 ) simply asserts that is a timelike one form of unit norm, as is usually true of 
a velocity held. The second of (1.6) asserts that the normal component of is not 
free but is determined by the other components of and the angle the membrane 
normal vector makes with a unit vector in the S direction. The third of (1.6) asserts 
that the velocity held has no component in the special S direction. The velocity one 
form held is specihed by p + 3 functions that live on the membrane. These p + 3 
equations are constrained by by the 3 equations (1.6) so that the data in the velocity 
is ehectively p functions of p + 2 membrane coordinates. 

The shape of the membrane and the velocity held described above are the vari¬ 
ables of our auxiliary dynamical system. As we will demonstrate below, the equations 
of this system turn out to be 

C/iXVi+ns(4-l)/S = 0 
P.‘(C'±V%) = 0 

where 

U^ = U-{U.n)n, U = dS+ nl{dS-u^dx^) (1.8) 

is the extrinsic curvature of the surface, the symbol dot (.) denotes contraction 
of indices (all indices are raised and lowered by the hat metric (1.3) ) and is 
the projector orthogonal to the 3 dimensional subspace spanned by the three one 
forms n^ and dS (so that projects onto a p dimensional subspace). (1.7) are 
a set of p -|- 1 equations that dehne an initial value problem for p -|- 1 variables (the 
shape of the surface and the constrained one form held n^). These equations are 
of hrst order in time derivatives of velocity but second order in time derivatives of 
the shape function, and so the data that determines a solution of these equations is 
initial shape of the membrane and its hrst time derivative together with the value of 
the velocity held at any instant of time. 

The equations (1.7) are the main result of this paper. These equations are com¬ 
pletely new and capture the nonlinear dynamics of the slow quasi normal modes of 
EST and collaborators (which previously had only been understood at the linearized 
level). We have learnt that the authors of the upcoming paper [18] have indepen¬ 
dently obtained a set of nonlinear equations that describe static large d membranes. 
While the membrane equations (1.7) of this paper apply to aribitrary dynamical con- 
hgurations, it should be possible to specialize these equations and thereby rederive 
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the results of [18] (at leading order in l/d) as a special case; we leave this exercise 
to future work. 

It has, of course, been previously noticed that the dynamics of black hole horizons 
resemble the equations of hydrodynamics on a membrane. The fact that the dynamics 


of black hole horizons is governed at large d by the precise equations (1.7) can be 


regarded as a mathematical realization of this idea in an appropriate context. In 
the title of this paper we have used the phrase ‘membrane paradigm’ to summarize 
our conclusion that black hole dynamics is precisely captured by an effective theory 
on a membrane. It would be interesting to investigate the relationship between the 
equations (1.7) and the work initiated in [19, 20] (see also [21] and subsequent work) 
that is more conventionally referred to as the membrane paradigm. 

As a test of our equations we have computed the shape function, normal vector 
helds, extrinsic curvatures and effective velocity one form helds for Schwarzschild 
black holes as well as Myers-Perry black holes (with rotations in upto [^^] planes) 
at large D, and have explicitly verihed that they obey the equations (1.7). We have 
also linearized (1.7) about the Schwarzschild solutions and computed the spectrum 
of small fluctuations; our results perfectly match the spectrum of quasinormal modes 
computed in [8]. ^ 

The reader might wonder how it is consistent to have a membrane fluctuating 
around in a theory of gravity without emitting gravitational radiation. We believe 
that the answer to this question is that our membrane does emit radiation, but 
that the amplitude of the radiation is of order any finite distance away from the 
membrane. The radiation may presumably be regarded as resulting from an effective 
coupling of the form 



Here h is the metric fluctuation to the exterior of the membrane. is the effective 
stress tensor of the membrane which is delta function localized on the membrane 
surface and must be conserved as a consequence of (1.7). It would be very interesting 
to derive this structure (and in particular the formula for in detail. We postpone 
this task to future work. 

The discussion of the previous paragraph would appear to suggest that our mem¬ 
brane loses very little energy to radiation. While this is possible it also possible that 
the smallness of the radiation held away from the membrane is actually a conse¬ 
quence of the extremely rapid fall off of a hnite energy radiation held due to the 
large d version of the the inverse square law. Restated, the highly dissipative nature 
of our membrane dynamics may be entirely due to “radiation” into the horizon (as in 

®Note that the quasinormal frequencies have imaginary pieces at the leading nontrivial order in 
the ^ expansion. It follows that the equation (1.7) are inherently dissipative, in contrast to the 
equations of hydrodynamics which are non dissipative at leading (perfect fluid) order. 
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the case of the fluid-gravity correspondence) but may also be partly due to radiation 
to inhnity. We postpone further discussion of this important issue to future work. 

The equations (1.7) are the principal results of this paper. Our equations sug¬ 
gest many interesting questions and generalizations. If the loss of energy to inhnity 
is small, it may be possible to recast (1.7) as the equations of conservation of the 
membrane stress tensor; this would yield a formulation of membrane dynamics anal¬ 
ogous to the usual formulation of hydrodynamics. It should be possible to imitate 
the analysis of [ 22 ] to uplift the area form on the event horizon to an entropy current 
on the membrane whose divergence is point wise positive for every consistent solution 
to the equations of motion. It would be interesting to hnd the hrst correction in ^ to 
(1.7) and also to hnd the analogue of (1.7) in the Einstein-Maxwell system. Finally 
it would be interesting to study the phenomenology of these equations; for instance 
to use them to study the head on collision of two spherical black holes. We leave all 
these questions to future work. 

2 Membrane like solutions in the large d limit 

2.1 Dimensional reduction of SO{d+ 1) invariant metrics 

As we have explained in the introduction, in this paper we study the vacuum Einstein 
equations (1.1), restricting attention to solutions that preserve SO{d+ 1) invariance, 
i.e solutions of the form (1.4). The ehective action for the metric and the scalar 
held (j) in (1.4) takes the form 

S = j [r + d{d - l)e-^ + ^^^^(<90)2^ (2.1) 

The equations of motion that follow from the variation of this action may be shown 
to take the form (see Appendix A) 

d d (2.2) 

2.2 A large d limit 

We are interested in ( 2 . 2 ) in the large d limit. A glance at ( 2 . 2 ) shows that 0 is 
not on the same footing as in this limit; derivatives of 0 are weighted with an 
additional factor of d compared to derivatives of the metric ® ^ It follows that 

^Recall that dV has a term involving one derivative and one Christoffel symbol, and so 
should be thought of as having one 0 derivative and one metric derivative at leading order. 

^The intuitive reason for this difference is clear; 0 controls the size of a d sphere. In the large d 
limit we should thus expect derivatives of 0 to be more ‘expensive’ than those of the metric. 
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nontrivial solntions in the large d limit must involve two length scales; a short length 
scale of order ^ for the variations of the metric and a longer length scale of order 
unity for the variations of (p. ® It is thus possible - and very natural - to choose 
coordinates for the equations in which the metric and 0 are both of order unity, but 
derivatives of 0 are of order unity while those of the metric are of order d. We choose 
to employ such a coordinate system in what follows. 

How can we describe a metric which varies on the length scale ^ in the large d 
limit, and yet is nontrivial over macroscopic length scales? We employ the following 
strategy. Choose a particular point Xq and blow up a region of size ^ about this 
point to unit size and then appropriately rescale the metric and dilaton gradient. In 
equations 

= 4 + 

= cPgah Hiii, = 'fo-lolgai = 

Xa = Va0 X d = a^v^0 


® We must then repeat the procedure of this paragraph for many Xq (so that the 
union of the patches about each of these points is the full manifold) and make sure 
that descriptions in distinct patches agree where they overlap. 

Focussing on the patch around Xq, (2.2) may be rewritten as 


1 

2 


VaX 


a 


= e-'^ 


d — 1 
d 



Rab a,Xb 


(2.4) 


where (2.4) is written regarding Gab as the metric (curvatures, Christoffel symbols 
and covariant derivatives are also constructed out of Gab-) 

At leading order (2.4) reduces to 


- 2G“'T;jXo = 4e-* - y 

2B„i = 


where 0 and Xa are constants. Note that d has disappeared from the leading order 
equations (2.5). 

gives a covariant estimate for the longer length scale; the square of the Riemann tensor 

gives a covariant estimate for the inverse fourth power of the shorter length scale. 

®We rescale V0 in order to account for the extra factor of d in derivatives of 0 as compared to 
derivatives of the metric. 

^°Note that 
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2.3 Flat solutions of the leading large d equations 

It is obvious that constant Gab is a solution to (2.5) provided 


= 4 ( 2 . 6 ) 

(here = G^^^XaXh) (2.6) is equivalent to (see (2.3)) 

e^{g^^d,ct)dA) = 4 (2.7) 

Its easy to check that (2.7) is satished at every point in (global) flat space. In the 
coordinate system of (1.3) this follows upon noting that 0 = 2 In S' and = 3“^. 


2.4 The black brane 

In this subsection we present a less trivial class of exact solutions to the leading 
order large d equations of motion (2.5). The solutions we present in this subsection 
maintain translational invariance in p + 2 out of the p + 3 dimensions in which we 
work; we refer to them as black brane solutions (in analogy with the black branes of 
asymptotically AdS Einstein gravity with a negative cosmological constant). 

Let Gf denote any metric whose components, {GF)ab are all constants. Consider 
the metric 

Gf + e-^0^ ( 2 . 8 ) 

where i? is a coordinate {dR is the direction in which translational invariance is 
broken) and O is a constant one form (i.e. each component of is a constant). It 
turns out that (2.8) is an exact solution to the equation of motion (2.5) provided 
that 

e‘^X.X = 4, (2di?-x).di? = 0, (2di? - x).0 = 0, 0.0 = 0 (2.9) 

(all dot products in equation (2.9) are evaluated using the metric Gf)- In Ap¬ 
pendix A.2 we have demonstrated that it is possible to choose coordinates in which 
the metric (2.8) takes the form 

= — 5 -, X = 2dR -\ - dX 

Xq Xq 

ds^ = 2dRdV - a (1 - dV^ + dY^dY^ + - 

At fixed xo it follows that (2.8) and (2.9) define a one parameter set of metrics 
upto coordinate transformations. In Appendix A.2 we demonstrate that the one 
parameter, a in (2.10), is a consequence of the scaling symmetry of the vacuum 
Einstein equations without a cosmological constant. 

^^Note that the metric (2.8) does not provide a completely unambiguous definition of the coor¬ 
dinate R and the one form O. The redefinition R ^ R + m and O —>■ O leaves (2.8) unchanged. 
It follows that the one form O and the coordinate R are defined only upto this ambiguity. 
^^Namely that uniform scaling of the metric by a constant maps solutions to solutions. 


dX‘^ 


— ax, 


0 


( 2 . 10 ) 
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The black brane solutions presented in this subsection are not really new. They 
can be obtained from the well known Schwarzschild solution (1.2) as follows. We 
choose any point on the horizon of the black hole, and then zoom into a patch of size 
1 of that point by performing the coordinate changes (2.3). The resultant metric 
turns out to be the black brane; the parameter a of the black brane is a function of 
tq of the original black hole. In other words the black brane solutions of this section 
are simply large d limits of appropriately ‘boosted’ versions of the ‘near horizon’ 
solutions of Schwarzschild black holes (of various sizes). 


2.5 Sewing black branes into a membrane 

In this snbsection we will present the construction of a class of metrics with several 
interesting properties. We start with flat space and perturb the metric in a manner 
controlled by a scalar function B and a null vector held O. 

Let B be an SO{d+l) invariant scalar function that lives on the hat space (1.3). 
SO{d + 1) invariance ensures that i? is a function only of the coordinates {w°‘,S) 
(see around (1.3) for dehnitions). We choose the function B so that its zeroes form 
a codimension one closed surface -for instance of the topology of a sphere times time 
- in the hat spacetime (1.3). We refer to this surface as the membrane. We require 
that B have only simple zeroes. The zeroes of B divide hat space up into two disjoint 
regions; we refer to the part of spacetime that includes inhnity as the outside, and the 
other region as the inside. We reqnire that B is negative on the inside and positive 
on the outside, and also that dB.dS > 0. 

In terms of B we dehne an auxiliary scalar function ^|J by the equation 


^/J 


1 I dB.dcf) 
2dB.dB 


( 2 . 11 ) 


It follows from the conditions on B that ij} <1 inside, that = 1 on the membrane 
and that -^ > 1 on the outside. 

In addition to B, let us consider a null SO{d -|- 1) invariant one form held O 
which lives in hat space. As O is null, 0.0 = 0 everywhere. We also demand that 


0 .( 


d(f) 

Y 


#)|£I.O = 0- 


( 2 . 12 ) 


where (2.12) is required to hold only on the membrane. The functions B and O, 
written as functions of (tc“, S) and (dtc“, dS), are assumed to be independent of d. 
Now consider the spacetime 




d'S'jiat + 


Ofj_0^dx^dx'^ 

rjpd+p 


(2.13) 


Here dsj-i^^ is the metric of hat space (for instance ds'ji^^ is given by (1.3) in a useful 
coordinate system). 
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The spacetime (2.13) has several interesting properties. First, it approaches 
the metric of flat space exponentially rapidly in the ontside region. The deviation 
of (2.13) from flat space scales like It follows, in particular, that (2.13) 

solves Einstein’s equations at large d at any point outside the membrane such that 
limrf^oo(V’ - l)d = oo. 

The only points outside the membrane at which it is not obvious that (2.13) 
solves Einstein’s equations is points for which -0 — 1 is scaled to be of order ^ as 
d is taken to inhnity. Recall that "0 = 1 dehnes the membrane surface. As 0 is a 
smooth function in the neighbourhood of the membrane, all such points lie within 
a distance ^ of the membrane. We refer to this neighbourhood of the membrane as 
the ‘membrane region’. We will now argue that (2.13) solves the leading order large 
d Einstein equations everywhere within the membrane region. 

In the membrane region ~ e(p+rf)(h-i)_ As 0 is a smooth regular d inde¬ 
pendent function in the neighbourhood of the membrane, it follows that 
and so the metric (2.13), varies over a length scale of order i. In contrast (p = 2lrLS 
varies over the length scale unity. 

Let us now study the metric in the membrane region in more detail, using the 
strategy outlined around (2.3). For this purpose it is useful to use 0 — 1 as one 
of our original coordinates, and choose R = d{'i/j — 1) a.s the corresponding rescaled 
coordinate for the blown up analysis of (2.3). Choosing in (2.3) to be any point 
on the membrane, it is easily seen that (2.13) reduces precisely to the black brane 
metric (2.8) at leading order in the large d limit. As the flat space metric and O 
were smooth functions in the original coordinates, they are constants after blowing 
up (to leading order in i) exactly as for (2.8). Moreover all the conditions (2.9) are 
also easily seen to be obeyed. The hrst of (2.9) follows from the fact that 0 = 2^5*. 
The second equation follows upon using (2.11) and the fact that B vanishes on the 
membrane. The third of (2.9) follows directly from (2.12) and the fourth of (2.9) 
follows form the fact that O is a null held. 

Finally, in Appendix A.3 we demonstrate that the event horizon of the metric 
(2.13) lies within the membrane region. For this reason we do not care whether 
(2.13) obeys Einstein’s equations inside the membrane (i.e. for 1 — 0 ^ i) as this 
region of spacetime is causally disconnected from the outside. 

In summary the metric (2.13) obeys the leading large d Einstein equations every¬ 
where outside its event horizon. This fact allows us to use (2.13) as the starting point 

for the construction of solutions to Einstein’s equations in a systematic expansion in 
1 

d- 

2.6 Equivalent starting points for perturbation theory and a ‘gange choice’ 

Recall that the metric (2.13) obeys Einstein’s equations only at leading order in large 
d. If two metrics of the form (2.13) coincide at leading order in large d (but differ at 
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subleading order) they should be regarded as equivalent starting points for pertur¬ 
bation theory. With this dehnition of equivalence classes we show in Appendix A.4 
that inequivalent classes of metrics (2.13) are labelled by the shape of the membrane 
and the value of the one form held O on the membrane. In other words functions B 
and one forms O related by 


0^0 + 60, B^aB 


(2.14) 


are equivalent starting points for perturbation theory. 

In order to proceed with our perturbative construction of metrics in the large d 
limit, we will hnd it useful to choose one representative in each equivalence class; we 
do this by imposing the following ‘gauge’ conditions on B and O 


{dB.V)V^B = 0 
dB.VO = 0 


(2.15) 


The hrst condition in equation (2.15) asserts that obeys the geodesic equation 
with affine parametrization. The second condition in (2.15) asserts that O is simply 
parallel transported along these geodesics. 

It is useful to note a simple consequence of (2.15) 

V^id^Bd^B) = 2d^BV^{V^B) = 0. (2.16) 

In other words the norm of the one form dB is a constant - not just along the 
geodesics, but everywhere in space. An overall scaling of B can then be used to set 
the norm of this one form to unity. With this choice of scaling, dB is, in particular, 
the unit normal one form to the membrane. We employ this choice in what follows. 
In the rest of this paper we use the notation 

?7,^(x) = df,B (2.17) 

which highlights the fact that d^B is the unit normal vector on the membrane. In 
this new notation (2.15) become 


n.Vn = n.VO = 0 (2-18) 

Finally let us turn to a more explicit parametrization of the one form O. Let 

0 = e\l,-Ua) (2.19) 

We assume SO vanishes on the membrane while a is nonzero but finite on the membrane. 

^^With these conventions there is an elegant geometrical construction of the function B in the 
neighbourhood of the membrane. Let us foliate spacetime with geodesics that pierce the membrane 
in a normal direction. The value of B at any point P is simply equal to the signed length along the 
appropriate geodesic from the membrane to P. 
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(where the hrst component in (1, —Ua) is the coefficient of dS, and the remaining 
components refer to the components of dw°'. It follows from the fact that O is null 
that 

= -1 ( 2 . 20 ) 

Moreover, the second of (2.15) implies that 

n.Vua = n.Vh = 0 (2-21) 

Finally recall that 0.{dip — ^) = 0 on the membrane. Using 
membrane, this condition reduces on the membrane to 

u.dB = dsB-^ (2.22) 

OsO 

Now it is easily verihed that the V^B'Vfj, derivative of the equation (2.22) vanishes 
(this follows from the second of (2.15)) implying that (2.22) holds everywhere, and 
not just on the membrane. 

3 Correction at first order in ^ 

In this section we will improve the metric (2.13) to ensure that it solves Einstein’s 
equations at hrst sub leading order in Our strategy to do this is straightforward. 
We will expand the metric (2.13) about Xq in the patch coordinates (2.3), but will 
keep both the leading and hrst subleading terms in the ^ expansion. In an appropriate 
choice of patch coordinates the metric reduces to the black brane solution (2.10) 
corrected by terms of order ^ that are linear in 6x^ = x^ — Xq. This metric fails to 
solve Einstein’s equations at hrst sub leading order in i. In order to hx this defect 
we allow the metric (2.13) to be corrected at hrst order in ^ and determine the form 
of the correction that ensures that Einstein’s equations are obeyed at hrst order in 
i. As we will see below, the correction is nonsingular if and only if the membrane 
shape and held obey the equations of motion (1.7). 

3.1 Interpretation of Taylor coefficients and constraints 

The functions B{x), h{x) and m^(x) can be Taylor expanded about any point Xq as 
follows 

Sx^Sx^ 

B{x) = B{xo) + V^B{xo)6x^ + V,V,B{xo)^^ + ... 

h{x) = h{xo)+ Vfj,h{xo)6x^ + ... 
u^{x) = u^{xo) + VuU^{xo)6x'' + ... 

where 6x^ = x^ — Xq. Let us define the functions 

V^5(x) = n^(x), Vf.'VuB^x) = K^,uix), V ^h{x) = hf,{x), V„u^{x) = u^y{x) 

(3.2) 


13 





It follows that the Taylor expansion (3.3) may be rewritten as 

B{x) = B{xo) + n^{xo)6x^ + — - -h • • • 

h{x) = h{Xo) + hf,{Xo)5x^^ + ... 

Uf,{x) = u^{xo) + u^u{xo)5x'' + ... 
n^{x) = n^{xo) + K^y{xo)5x'' + ... 

Recall that the norm of n^{x) is a constant (see (2.16)) and that we have chosen to 
normalize B to set this constant to nnity. It follows from the hrst of (2.15) that 

n^^{x)K^,{x) = 0 (3.4) 

Using (3.3) and (3.4), it follows that the extrinsic cnrvatnre at any point Xq on the 
membrane snrface is given by simply Kf^i^{xo). 

The derivatives of the velocity held are not all independent; differentiating (2.22) 
we hnd 

(3-5) 

and differentiating (2.20) we have 

= 0 (3.6) 


3.2 Data 

At any point on the membrane, there is a p dimensional snbspace orthogonal to the 
three dimensional snbspace spanned by In this brief snbsection we list a 

basis of independent Taylor expansion coefficients, labelled by their transformation 
properties nnder the local tangent space SO{p) gronp. Recall that the Taylor ex¬ 
pansion coefficients u^y and are not independent as they are constrained by 
(3.4), (3.6),(2.21) and (3.5); these eqnations can be nsed to determine all the ‘other’ 
expansion coefficients in terms of the independent coefficients listed below 

• Scalars: 

= Kss, S2 = S3 = u^u'^K^y, = V^'^K^y 

S3 = K, S, = S-r = V^''u,y 

• Vectors: 

= V^^u’^u^y, u/ = Vi^ = 


Symmetric Tensor: 


_ 'pfiap^p 


K, 


a0 


Va0 
P 


9(t> 


_ p^^apfJ,j3 


2 p 
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• Anti-symmetric tensor: 


_ 'pfia pul3 '^a/3 Ufjg 


As we have explained above, onr ‘gange’ condition (2.15) ensnres that normal 
derivatives of h and the velocity held vanish (see (2.21)). This ensnres that the 
derivatives of these qnantities are eqnal to derivatives projected to the world volnme 
of the membrane, i.e. 


M.V = M.V — (n.M)n.V, ds = ds — ns{n.'V) 


(3.7) 


As the h and velocity helds are physical only on the snrface of the membrane, all 
derivatives in the data above (and in the hnal hrst order metric for onr black hole 


below) shonld actnally be thonght of as projected derivatives. With this replacement 
onr hnal resnlts for the metric and eqnation of motion apply withont reference to 


the arbitrary gange choice (2.15). 

3.3 Structure of the first correction to the metric 

As mentioned above, we will need to correct the metric (2.13) in order to ensnre that 
Einstein’s eqnations are obeyed at hrst snbleading order in i. We adopt S, the size 
of the d sphere as one of onr coordinates. This choice leaves ns with p + 2 remaining 
coordinate freedoms. We hx this freedom in the following manner. Let dijj, dS — Sdip, 
O and Vj^dx^" form a basis of the p + 3 dimensional space of one forms at any point. 
The correction metric can be regarded, point wise, as a symmetric qnadratic form 
of these one forms. Onr choice of gange is that the only terms involving d^|J in the 
correction metric are of the form 


d'ip(^{d-il: - —); 


The correction metric will not contain terms proportional to dijj^d'ijj, d'ljj^O and 
dip (8) Vjpdx'^. 


With this choice of gange the metric (2.13), corrected to hrst order in i, mnst 
take the form 


ds^ = p^^dx^^dx'' + ip-^'^+P^O'^ 


+ ^ [iLi(a:“)02 + 2K2{x^){dS - SdiP) ® O 
d L 

+ K^{x^){dS - SdiPf + 2K^{x’^)diP ® {dS - SdiP) (3.8) 

+ 2Ql{x^)0;P^dx>^dx'^ + 2Qf (x“)pfda:^(d^ - SdiP) 



use the symbol ® to denote the symmetrized product of two one forms. If A = and 

B = B^dx^ then A® B = A^B^dx>^dx‘'. In our convention A^ = A ® A. 
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The O(^) correction piece in (3.8) has been added to ensure that (3.8) solves 
Einstein’s equations at hrst subleading order in This correction is needed because 
the hrst line of (3.8) does not, by itself, have this property. Let us recall why this is 
the case. Upon expanding the hrst line of (3.8) in a patch to hrst subleading order 
in i, we hnd the black brane background corrected by 0{l/d) huctuations. Now the 
huctuations - which are all proportional to the data of the previous subsection plus 
a constant piece fail to solve the leading order Einstein equations, giving rise to 
sources proportional to hrst order data. Moreover the black brane metric itself is a 
solution only to the leading order Einstein equations, and so yields a constant source 
at order ^ when plugged into the exact Einstein equations. As the correction part 
of the metric in (3.8) is chosen to cancel these sources, it follows that the unknown 
functions in (3.8) must take the form 


Ka{x^) 




where R 


Y,K^{R)Sa^ a = {1,2,3,4}, = 1 

^=0 

« = {i.2} 

A=1 

2 

A=l 

d{i/j — 1) 


(3.9) 


The unknown functions Q^{R), Ta{R) are now functions only of R and 

can be regarded as constants in the other black brane directions within a patch (their 
variation in these directions ahects the metric only at 0{l/d‘^)). 


3.4 Results 

In order to determine the unknown functions in (3.9) we now expand the metric (3.8) 
in patch centred about a particular point Xq and obtain the patch metric expanded to 
hrst order in ^ (including terms that come from the Taylor expansion of the hrst line). 
We then move to a convenient set of patch coordinates and plug the resultant metric 
into Einstein’s equations. As we explain in detail in Appendix B, we hnd diherential 
equations for all unknown functions in (3.9) which we proceed to solve. The solution 
to these equations turns out to include a singularity within the membrane region 
(see Appendix B) unless the equations (1.7) are obeyed. When those equations are 

^®The constant arises from the fact that reduces to e~^ only at leading order at large d; 

there are corrections to this formula at first subleading order in i. 

^^The appearance of S'o = 1 is a consequence of the fact that we have constant sources (not 
proportional to the first order data of the previous subsection) as explained above. 
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obeyed the unknown functions in (3.9) are all quite simple and we find that the 
hnal metric ((3.8) with explicit solutions for the functions in (3.9)) is given by 


= rjujydx^dx’^ + -0 ^^^'^\0^dx^Y 
^-d 


{0^,dx^) 


d ^ ^ 
2 


ibi(a:")(0^da;^) + 2K2{x°){dS - Sdi^) - Kv{x^){V^up,)dx'' 


O 


where 


Ki{x^) = - 


-c® + llc^ - 23 c2 + 11 

(c2 - 2f (c2 - 1) 


^ {c^ — 1 ) (3c^ — 8 ) 


R- 


c2i?2 + 2 (-4c^ + 7 c 2 + 1) i?' 
2 c 2 (2 - . 

+ (-4c^ + 6 c 2 + 2 ) R: 


2 (c 2 - 2Y 

{u^u^K^pYSns) 


R^ 


C2 (1 - C2) (2 - C2) 


(uf^KsYiSns) 


+ 


R 


1 — c" 


125%/! - (r'“'K^MSns)] 


(3,10) 


(3.11) 


^®In this paper we have performed all our computations only for the special case p = 2. The 
formulae presented in equations (3.11), (3.12) and (3.13) are accurate only for p = 2. We do 
not, however, need to perform a new computation in order to obtain the analogues of these three 
functions at arbitrary p. This is because a solution at p = 2 may be viewed special case of a 
solution of the membrane equations at all larger values of p - simply by moving some of the zm 
coordinates into the set of Wa coordinates. This operation changes the meaning of S', and so induces 
a non-trivial redefinition of the velocity one form (recall u is defined so that O oc (dS — u)). This 
requirement of ‘redistribution invariance’ together with the explicit p = 2 expressions presented 
below in (3.11), (3.12) and (3.13) is sufficient to completely fix the three functions Ki, K 2 and 
at all p. We postpone the presentation of the explicit expression for the first order correction to 
the membrane metric at arbitrary values of p to an upcoming paper. Similar logic may be used to 
constrain the p dependence of the membrane equations of motion (1.7). In this case we find that 
the requirement of redistribution invariance infact forces the membrane equations to take the form 
(1.7) at all p. This fact will also be explained more fully the upcoming paper mentioned above. 
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h\ - 6c^ + 1) + (c® - 3c^ + 2c^) R 
) = ® ^ ^-Vti—7T-^ 




(3c^ — 6c^ + 1) + (c^ — 2) c^R 

[ c4 (2 - c^f . 



(3.12) 


2(2 - c^)idsh) - 2(1 - ^){u^deh) 

- {V^^K>.u){ns) + (1 - c2)(P^X,) 



(3.13) 


Here 


= The projector perpendicular to and dS 


2 1 2 
c = 1 — n _5 

R = d{'il) — 1) 


(3.14) 


All raising and lowering of the indices have been done using the flat metric rj^y. 

(3.10) is our hnal result for the metric of membrane spacetimes corrected to first 
order in 2. We recall that (3.10) solves Einstein’s equations to the relevant order if 
and only if (1.7) are obeyed. 

4 Stationary black hole solutions 

In this section we provide some examples of exact solutions to Einstein’s equations 
that reduce to membrane spacetimes at large d. We explicitly verify that these 
solutions obey the equations (1.7). 

4.1 The Schwarzschild black hole 

The metric for the Schwarzschild black hole of unit radius is given by 


ds'^ = —dt^ + dS"^ + dxidx^ + 



Hence the zero norm one form O is given by 


u^dx^ 


O 



(4.2) 


18 






















(4.3) 


where = XiX^. The membrane surface is dehned by 

S'^ + XiX^ = 1 


the normal to the surface which has unit magnitude on the surface is given by 


SdS + rdr 
n = — , 

\/52 + r2 

2 

We see that the conditions u.u = —1 and u.n = n® — ^ are trivially satished. 

Let us now focus on the membrane surface S'^+r^ = 1. The non-zero components 
of extrinsic curvature of this surface are obtained as 



Kss = l- 
Ks^ = -Sxi = K,s 

Kij = 6ij - XiXj (4.5) 

We will now check that the Schwarzschild black hole obeys the equations (1.7). 
The value of the constant in these equations turns out to be when evaluated 
on this solution. The various terms appearing in the scalar equations evaluate to 

{2-cyKss = = ^V, 2{l-c^){2-c^)Ksu = 2il-S^)S\^ 

(4,6) 

Plugging all this in the scalar equations of (1.7) we see that it is satished. The two 
terms in the vector equations of (1.7) individually evaluate to zero on the membrane 
and hence it is also satished. 

In conclusion, the Schwarzschild black hole spacetime takes the membrane form 
at large d; the shape of the corresponding membrane and its velocity held solve the 
membrane equations of motion. 

4.2 The Myers-Perry black hole 

Consider the Myers-Perry black hole [23, 24] with independent rotations in q distinct 
two planes. In order to describe this solution within our formalism we need to choose 
p -|- 1 > 2g. We hnd it convenient to saturate this inequality so that p -|- 1 = 2g or 

In Appendix A.5 we demonstrate that the Myers-Perry black hole with rotations 
in 272. planes is described by the metric 

P+1 

= —dt^ + ^(dXi -I- dy^) + dS^ -\ ^ (4.7) 

i=i P 
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where = zmz^ as around (1.3) and 


P+1 


p+i 

2 


dp-Yl 


i=l 




{xidyi - Vidxi 


( 


A = 


\ 


m 


P+l „2 / P +1 „ 2,,2 

n.ii(i + ^) (1-E.il ^ 




(4.8) 


and 


E 4 ^ + ^ = 1 - 

^ + of p 2 

.=1 ‘ 

2 _ + Si 

P^ + a?' 


(4.9) 


SO that 


' P+1 
2 


pdp = 


Z^7 = l rn2 4-27^'l2 W „4 


E- 


Xidxi + Uidui , S' 


2=1 


p" + < 




(4.10) 


^i-1 (p^+a|)2 ' p-* 

The metric (4.7) takes the form (2.13) with = p (recall p is dehned implicitly by 
(4.9)). It is not difficult to verify that Ofj,0^ = 0 on the membrane p = 1. It is not 
difficult to show that 


O (X dS — u^dx^ 


where 


p+1 

2 


„ = I dt + ^ 


2=1 


Xi H" diUi j Hi diXi 

dXi -f- —^^— dy-i 


1 + af 


1 + af 


« = i-E 


p+1 


2=1 


1 + of 


(4.11) 


It is easily verihed that = —1. The unit one form normal to the membrane is 
proportional to dp; upon normalizing we hnd 


n = 


p+1 

E.ii 


( 1 +^ 


^2 


' p+1 

E 

2=1 


Xidx"^ + yidy^ 
1 + af 


+ SdS 


We have explicitly checked that 

n.n 

u.n = ns - 

ns 


(4.12) 
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so that (2.22) is obeyed. It follows that the metric (4.7) agrees with our form for the 
membrane metric at leading order in large d. 

We have also verihed that the expressions (4.11) and (4.12) obey the equations of 
motion (1.7). In order to verify the scalar equation (hrst of (1.7)) we explicitly com¬ 
puted the extrinsic curvature tensor. Using our explicit (and lengthy) expressions, 
we hnd that the scalar equation (hrst of (1.7)) is indeed satished. 

We have also verihed that the vector equation (the second of (1.7) is obeyed. 
Let us specialize to the case g = 1 or p = 1. On the membrane surface we hnd 


Va {(ns) ~ (n.n)?7,.V) Ub) = 


= Vl ((1 + n|) {ds - nsn.V) Ub) (4.13) 

(where we have listed the S, t, x and y components of the vector in that order). It 
follows that the second of (1.7) is satished. Note that each individual term in this 
equations is non-zero; however the two terms cancel between each other. Although 
we do not explicitly report our results here, we have also verihed that the vector 
equation is satished for the case g = 2 or p = 3. 

5 Linearization of membrane equations about the Schwarzschild 
solution 

In the previous section we verihed that rotating black holes at large d are solutions 
to our membrane equations. The equations (1.7) apply to more than stationary 
conhgurations; they completely capture time dependent black hole dynamics as well. 

In this section we apply our equations to a non stationary problem. We determine the 
most general solution to the equations (1.7) in linearization about the Schwarzschild 
black hole solution. The solutions we hnd should map to the most general regular 
linearized solution of general relativity (with frequencies of order unity ) about the 
Schwarzschild black hole solution, i.e. the ‘light’ quasinormal modes determined by 
EST (see the introduction). Indeed we hnd that the spectrum of the solutions we 
obtain agrees precisely with the spectrum of the light quasinormal modes of EST to 
appropriate order in 4. We regard this match as a signihcant check on our membrane 
equations (1.7) 

In the previous section we demonstrated that the membrane surface 

S^ + r^ = l (5.1) 


/ 0 \ 

a2(a2+l)(x2+y2^(^2a'‘-2a2(a;2+y2_2)-a;2-y2+2) 

(a^—a2(x2+y2_2)+l)^ 

a(a^+l)'^j;(2a'*—2a^ (x^+y^— 2 ) —y^+2) 

(a^—a2(x2+y2_2)+l)^ 
a(a2+l)'^x(2a'*—2a2 (^x'^+y'^—2^ —x'^—y'^+2^ 

\ (a^—a2(a:2+j/2_2)+l)^ / 
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(recall w'^Wi where the index i runs only over spatial coordinates) together 

with the velocity one form held on the surface of the membrane 

n = df + (-^) dr (5.2) 

solves the membrane equations (1.7), and that this solution corresponds to a Schwarzschild 
black hole of unit radius. In order to study small huctuations about this solution we 
let the membrane huctuate, i.e we let our membrane surface be given by 

5-Vl-r2-a/(r,f,d0 = 0 (5.3) 


Equivalently, the S coordinate on the membrane is determined in terms of r, d, v by 
the equation 


We also allow the velocity held to huctuate, i.e we set 


(5.4) 


u = + a5u^ dt + dr + aduidO'^ (5.5) 


where S is given by (5.4). The parameter a in (5.3) and (5.5) measures the amplitude 
of the perturbation; through this section we work at linear order in a. 

Our membrane ehectively propagates in the p + 3 dimensional space with metric 


ds^ = —dt^ + dwidw^ + dS‘^ = —dt^ + dr^ + r‘^dVt^ + dS"^ (5.6) 

The normal to our perturbed surface is given, to 0(a), by 
n = (\/l — + ar(l — r'^)drf)dS — ay/l — r'^dtfdt + (r — a\/ (1 — r'^)^drf)dr 

—aVl — r'^difdO'' (5.7) 


Once again the variable S in (5.7) is given by (5.4). Here d* are an arbitrary set of 
coordinates on the unit p sphere in the metric (5.6), and we have normalized n to 
ensure that n.n = 1 upto corrections of order a^. 

2 

Recall that the velocity held u is constrained to obey u.u = — 1 and u.n = ns — ^. 
These conditions allow us to solve for 6ur and 6ut in terms of /; we hnd 

= ( 1^-n 2^ + ^9rf)) 

r(-l + r0 

Sut = rdur — af (5.8) 

It follows that we have p + 1 independent huctuation helds (/ and Ut). 
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5.1 The scalar equation 


We find it convenient to use t, r, 9^ as coordinates on the membrane. The dynamics 
of / is controlled by the scalar equation (the first of (1.7)). Evaluating the extrinsic 
curvatures associated with the normal vector field (5.7) (we performed the relevant 
computations on Mathematica), we find that / obeys a linear partial differential 
equation of second order in derivatives in t and r. The equation has no derivatives 
in the angular coordinates 0*. 

It is useful to expand the fluctuation function / as 

/(r, t, 9^) = (5.9) 

Im 

Here I and m refer to labels for scalar spherical harmonics of SO{p + 1). Hs a label 
for the representation of SO{p + 1) and m is a collective label for the quantum 
numbers of states within a given representation. 

Plugging this expansion into the linear equation for / yields an ordinary dif¬ 
ferential equation of Gim{r). Because the equation for / had no derivatives in the 
angular directions, all functions Gim obey the same equation. Making the replace¬ 
ment Gim —)■ G and ojim —)■ cu we find that the equation takes the form 

— -I- w{—2i -I- tc) -I- r^(l -|- 2iw — 2w'^))G{r) 

+r{ — l + r^)(2(r^ + iw — ir‘^w)G'{r) + r(—1 -|- r^)G"(r)) = 0 (5.10) 

where ' denotes derivative w.r.t. r. The two linearly independent solutions to the 
equation (5.10) are given by 

Gir) = ^ ^ 

-|- r‘^){i -f Aw) ^ ' 

Regularity of these solutions at r = 0 requires the power of r in (5.11) to be a non 
negative integer. One of the two solutions above is proportional to Y provided 

ws{j) = -i{j - 1) ± \/j - 1 (5.12) 

While j is required to be a non negative integer this condition is in itself not sufficient 
to ensure regularity for several different reasons. 

Let us first examine the lowest allowed value of j, namely j = 0. At this value 
of j the two frequencies (5.12) are tcs(O) = 0 and tcs(O) = 2i. It is easily verified 
from (5.8) that the mode with tCs(O) = 0 has 6ur = 0 and 6ut regular everywhere. 
Consequently this mode is a regular small fluctuation of our membrane equations. Its 

we measure highest weights of representations by eigenvalues under rotations of independent 
two planes, the representation labelled by I has highest weights (/, 0, 0 ... 0). 

Which one depends on our precise definition of the square root. 
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physical interpretation is simple; it corresponds to rescaling the black hole. When, 
on the other hand, tCs(O) = 2i, we hnd from (5.8) that 6ur oc It follows 

that the velocity held corresponding to this mode is irregnlar at r = 0 , so this mode 
is irregnlar and so nnacceptable. This is reassnring as w = 2i wonld have represented 
an instability, contradicting the well established stability of Schwarzschild black holes 
to small hnctnations in all dimensions. In snmmary, the only allowed freqnency at 
j = 0 is ws = 0. 

Let us next turn to j = 1. At this value of j the quadratic equation for ^^(l) 
yields a double zero; both roots in (5.12) are tc = 0. The fact that we hnd a double 
root indicates that j = 1 is a special case that has to be dealt with separately; the 
two allowed time dependences at this value of j are actually / constant in time and 
f (xt These two ‘zero’ modes have a simple physical interpretation. The constant 
mode corresponds to an inhnitesimal translation of the black hole, while the mode 
linear in time corresponds to an inhnitesimal boost for the black hole. 

The next issue concerns the angular dependence of modes. Recall that the 
angular variation of / did not enter our equation, so the reader may at hrst suspect 
that we have an inhnite degeneracy of modes at every j, corresponding to the inhnite 
number of possible angular variations, seemingly unconnected with radial behaviour. 
This is not the case. A scalar function of the form r^Yim{0) is regular at r = 0 only 
if 23 

j = I, / + 2, / + 4 .... 

2 ^ It follows that we hnd a regular solution for Gimir) provided 

w = ws{i), j =/, / + 2, /+ 4... (5.13) 

In particular the minimum frequency for Gim{r) is Ws{l)- 

In summary the spectrum of frequencies of oscillation in the scalar sector is given 
by 0Js{j) for j = 0,1, 2,3,.... For any given value of j, the modes that have frequency 
oJsij) are scalar spherical harmonics with 


I -2,j -4... 


In other words modes with diherent angular momenta share the same frequency. 
Degeneracies are unusual in physics; when they occur they are usually explained by 
a symmetry. In the current context the relevant symmetry is obvious; hnctnations 

is easily verified that for j > 1 Sur is well behaved at r = 0. 

^^Note that t = Ine*; the linear dependence in time follows from the application of a standard 
rule. Recall that we plug ansatz a;“ into a second order linear differential equation and find a double 
root for a at a = Oc then a basis for the solution space of the equation is given by and Inx. 

^^The underlying reason for this is that spherical polar coordinates are singular at r = 0. 

^^This is most easily seen by noting that an SO{p + 1) spherical harmonic of degree I is the 
restriction of a polynomial of degree I (in the p + 2 embedding Cartesian variables) to the unit 
sphere. 
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about the Schwarzschild black hole appear in representations of SO{p + d + 2). The 
different SO{p + 1) spherical harmonics that descend from a single SO{p + d + 2) 
spherical harmonic must have the same frequency. 

Now consider an SO{p + 1 + d + 1) spherical harmonic that transforms in the 
representation (j, 0, 0 ... 0) (as usual we label representations by the highest weights 
under rotations in the independent embedding two planes). This representation can 
be decomposed into a sum over irreducible representations of SO{p+ 1) x SO{d+ 1). 
The number of representations that appear in this decomposition is very large. If 
we are interested only in those representations that happen to be singlets under 
SO{d +1), it is not difficult to convince oneself that the branching rule is 

[§] 

(j, 0, 0 ... 0) - 2m, 0, 0 ... 0) (g) J 

m=0 

In other words the set of SO{p + 1) spherical harmonics that have frequencies ujj 
are precisely obtained from a single angular momentum j spherical harmonic of 
SO{p + d + 2) subject to the constraint that we restrict attention only to SO{d + 1) 
singlets. 

In summary, the spectrum of scalar fluctuations is given by 0Js{j) for non negative 
j. Eigenmodes with frequency ujs{j) are given by the projection onto SO{d + 1) 
singlets of an angular momentum j spherical harmonic of SO{d + p + 2). This result 
is in perfect agreement with the frequency of light scalar quasinormal modes about 
the Schwarzschild black hole obtained by EST in Equation (5.29) of [8]. 

5.2 The vector equation 

In the parametrisation of this subsection the vector equations of (1.7) become 

(-2 + r'^)ui + r(l - r‘^)drUi + (1 - r^)dtUi = -(1 - r^)dif (5.14) 

where Ui is a one form held that lives on the unit p sphere. 

(5.14) is an inhomogeneous or forced linear differential equation in the one form 
variable Uj. The solutions of any such equation take the form H + P where P is 
a particular solution and H is the general solution to the homogeneous part of the 
equation (i.e. (5.14) with the RHS set to zero). 

The particular solution may easily be seen to be regular (the key point here is 
that the frequencies of the homogeneous part of (5.14), which we work out below, 
never coincide with the scalar frequencies of the previous subsection) and is of no 
particular interest to us; we do not bother to work it out. Our interest is in the 
spectrum of the homogeneous solutions, to which we now turn. 

is natural to choose the particular solution of this equation as a linear combination of terms 
with the scalar frequencies listed in (5.12). At any given frequency, the particular (or forced) solution 
is proportional to the coefficient of the scalar mode in / at the same frequency, and multiplies a 
regular function of r. 
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The most general solution to the homogeneous equation 

(—2 + r^)ui + r(l — r‘^)drUi + (1 — r‘^)dtUi = 0 (5.15) 

is easily determined. It is given by a linear combination of terms of the form 

M^(r) = 

Here qi{9) is an arbitrary one form held on the unit sphere. 

The solutions above are clearly irregular at r = 0 unless 2 — is a non negative 
integer. This is a necessary and not a sufficient condition for regularity. 

In order to understand the sufficient conditions for regularity more clearly it is 
useful to decompose 

qi = Vi + diX, where VjU* = 0 (5.16) 

and A is a scalar function on the sphere. Vi may then may be decomposed in a basis of 
vector spherical harmonics, while A may be decomposed in a basis of scalar spherical 
harmonics. 

Let us dehne the vector spherical harmonic as the set of vector functions 
on the sphere that transform in the (/, 1, 0,... 0) representation of SO{p + 1) (we 
label representations by eigenvalues of the rotation operators in the independent two 
planes in the embedding p+1 dimensions). It is not difficult to convince oneself that a 
one form held whose angular part is proportional to the vector spherical harmonic 
is regular if and only if the function that multiplies it is a linear combination of terms 
of the form where 

J = /y + 1, /v + 3, /y + 5 . . . . 

Let us now turn to the regularity of dX. As we have already noted in the last 
subsection, a scalar function A is regular if and only if the coefficient of the scalar 
spherical harmonic (dehned as the set of functions in the representation (Z^, 0, 0 ..., 0) 
of S'0(p + 1)) is a linear sum of terms of the form where 

j = hJs + + 4: ... 

Notice that for /^ = 0 dA = 0 so that the minimum value for Is is unity. 

Putting all these facts together, we hnd that we have regular solutions to the 
homogeneous vector equation at 

^ = ^v{j) =-i{j j = l,2... (5.17) 

This spectrum is in perfect agreement with Equation (5.21) of [8]. 

can also list the eigen functions. The allowed values of j for the frequencies of the 
vector spherical harmonic are j = /y,Zy + 2,Zy+4,... while the allowed j for the frequencies 

for the gradient of the scalar spherical harmonic are j = Zg — 1, Zg + 1, + 3 . We several 

degeneracies in the frequencies at any particular value of j. As in the previous subsection we 
presume that all representations with a given wy(j) have their origin in a single vector spherical 
harmonic representation of SO{d + p + 2), but we have not carefully verified this claim. 
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Note that wy vanishes at j = 1. This mode corresponds to turning on an in- 
hnitesimal rotation of the black hole (recall that rotations, like angular momenta, 
transform in the representation (1,1,0...0) - i.e. the lowest vector spherical har¬ 
monic representation - or the appropriate orthogonal groups. All modes with j > 2 
have an imaginary part to their frequency that corresponds to an exponential decay 
in time. It follows that the Schwarzschild black hole at large D is stable to vector 
fluctuations (as it was to scalar fluctuations) in accordance with well known results. 

Acknowledgements 

We would like to thank P. Ajit, B. Craps, Y. Dandekar, S. Das, R. Emparan, A. 
Ghosh, R. Gopakumar V. Hubeney, K. Imbasekar, S. Mandal, M. Mandlik, S. Nad- 
karni, P. Nayak, M. Rangamani, A. Sen, N. Seiberg, S. Trivedi, S. Wadia and A. 
Yarom for useful discussions. The work of S.B. was supported by an India Israel 
(ISF/UGC) joint research grant. The work of A.D is supported by an Inspire Fel¬ 
lowship from DST, India. S.M. would like to acknowledge the hospitality of ICTS 
where this work was intitiated and of the Harish Chandra Research Insitute, USER 
Kolkata, the University of Swansea and the Galileo Galilei Institute while this work 
was in progress. The work of S.M. was supported by an India Israel (ISF/UGC) 
joint research grant. We would all also like to acknowledge our debt to the people of 
India for their steady and generous support to research in the basic sciences. 

A Technical details 

A.l Simplification of the dimensionally reduced equations 

Upon varying (2.1) w.r.t the scalar 0 we hnd the equation of motion 

R + e-^{d - 2)(d - 1) - ~ ^\ d0)^ -{d- 1)VV = 0 (A.l) 

Varying w.r.t. we hnd 

(A.2) 

(A.2) may be somewhat simplihed by subtracting half of (A.l) multiplied by and 
the equations of motion may more simply be written as 

R + e-^{d - 2)(d - 1) - -{d- 1)VV = 0 

^ (A.3) 

((<i - l)e-* - - Ivv) 
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By subtracting the trace of the second equation from the hrst we hnd equivalently 


e-*{d - 1) - J W)" - ivV + ((“i - !)<= 

Rflu = ({d- l)e"'^ 


((<i - l)e-* - - ivv) 

d - l)e-* - - ivy) 


0 


(A.4) 


where the indices /r, u run over p + 3 values. Notice that the two terms in the 
scalar equation are simply proportional to each other, so that the scalar equation is 
equivalent to 



Given this condition, the second term in the vector equation vanishes, and our system 
of equations reduces to 



(A.5) 


A.2 Details of black branes 

A.2.1 Standard coordinates for the black brane 

In order to verify that (2.8) is indeed a solution, it is useful to move to a convenient 
coordinate system. Let xq denote the (constant) radius of the d sphere about the 
point of interest so that 



We then dehne the coordinate X by the requirement 



Xq 

We also dehne the coordinate V by the requirement that 

dV = 0 

These dehnitions are sensible because the components of y and O are constants 
At this stage we have three special coordinates; R, X and V. Let the remaining 
p coordinates be denoted by y*. We employ coordinate redehnitions of the form 
^ + C\iR + CyV + CxX to remove all cross terms between y* and (A, V, R). 

We then employ redehnitions of the form y* = Qy^ to set the metric in the F* 

emphasized in the previous footnote, this dehnition fixes V only upto a scaling. We will hx 
the ambiguity of the scale factor below. 




directions to 6ij. With all these conventions we now proceed to display an explicit 
form of the metric ( 2 . 8 ). 

The choices we have made in the previous paragraph clearly ensure that the 
metric is a direct sum of dY^dV^ and a three dimensional metric involving dX, dV, 
and dR. We will now determine the form of this three dimensional metric. 

The conditions (2.9) ensure that there are no cross terms between dX and either 
dV or dR. It follows that the metric takes the form 

-+ bdR^ + 2CdRdV - DdV^ 

1 — axQ 

where a, b, C and D are arbitrary constants. The fact that dV is null sets b to 
zero. As we have emphasized in the footnotes above, the coordinates V and R are 
ambiguous upto a coordinated rescaling and shift. We choose to £x this ambiguity 
by setting C = 1. Making this choice, the condition on the squared norm of y sets 
D = a, and the metric (2.8) takes the explicit form 


e '^ = — 5 -, y = 2dR H- dX 

Xg Xq 

ds^ = 2dRdV - a (1 - dV^ + dY^dY^ + 


dX^ 

1 — axg 


(A. 6 ) 


We have used an additional shift in R to scale the factor a outside e ^dV^ to ensure 
that the event horizon lies at i? = 0 . 


A.2.2 Uniform scaling of black branes 


The black brane metric (A. 8 ) appears in a two parameter family labelled by Xq and 
the parameter a. There is a simple explanation for the parameter a; its existence 
follows from the fact that Einstein’s equations in p + d + 3 dimensions are invariant 
under an overall rescaling of the metric. In terms of the reduced fields described 
above, the map 

(fi'/i!/) e'^) (Agf^y, Ae'^) (A.7) 


takes solutions to solutions. It is easily verihed that this rescaling operation does 
indeed map black branes to black branes. More concretely, the rescaling operation 
(A.7) applied to the black brane (A. 8 ) yields a new black brane 


e ^ = 


X, 


2 ’ 


y = 2dR H- dX 

Xo 


= 2dRdV - d (1 - e"^) dV^ + G^jdY^dY^ + 




(A. 8 ) 


axn 


with 

d=^, xl = Axl, V = AV, X = \/AX (A.9) 

XT. 

The important point here is that this rescaling changes the value of a, and so gener¬ 
ates black branes with all values of a starting with black branes at a = 1 . 
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A.3 Location of the event horizon 


We will now demonstrate that the spacetime (2.13) has an event horizon located 
‘near’ to the surface '0 = 1. As usual we consider a solution of Einstein’s equations 
that eventually settles down into a rotating black hole solution. In this context the 
event horizon is the unique null manifold that reduces, at large time, to the event 
horizon of the eventual black hole. In order to compute the event horizon manifold 
of the spacetime (2.13), let us set up a coordinate system in the neighbourhood of 
the membrane. Let 0 be one of our coordinates and let dehne the remaining p + 2 
coordinates. Let the event horizon manifold be given by the equation 

'tp = i}H{z‘^) (A.10) 


We will now demonstrate that 




01 

l + ^ + 

a 


The argument is very simple. The event horizon is determined by the requirement 
that the one form 

dtp + datpHdz°‘ 

is null. Let us work in a small patch around the point Xq which lies on the horizon. 
We have argued above that the metric is that of the black brane with the role of tp 
played by the coordinate R. Ignoring derivatives the null surface lies at the location 
e~^ = a. While derivatives are nonzero, they are of order i. It follows that the 
true null surface deviates from e~^ = a only at order i. In other words we have not 
merely established our result; we have also determined the precise location of the 
event horizon at leading order in it is simply the event horizon of the local black 
brane. 


A.4 Equivalence in data off the membrane 

In this short appendix we argue that the replacement 

0^0 + 60, B ^aB 

where 60 vanishes on the membrane, and a is any function that does not blow up 
or vanish on the membrane, generates equivalent metrics of the form (2.13). 

In order to see this we note that in the large d limit 

IpP+d 

As tp moves away from unity, in other words, this function dies off extremely rapidly. 
In other words the metric (2.13) only cares about the vector held O in a neigh¬ 
bourhood of the membrane where 0 — 1 = 0{l/d). In this region the hrst Taylor 
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expansion of O - and the second Taylor expansion of "0 - away from the membrane 
each contribnte to the metric at order To leading order, in other words, the metric 
cares only abont the location of the membrane, O on the membrane, and dijj on the 
membrane. 

We will now argne that (i0 on the membrane is also completely determined by 
the shape of the membrane, dip is proportional to the the normal vector held of the 
membrane snrface. As a conseqnence it is hxed by the shape of the snrface npto 
a position dependent normalization. Now the normalization of ip is nnambignonsly 
determined by the reqnirement {Sdip — dS).dS = 0. Indeed onr dehnition of ip 
antomatically obeys this eqnation. As the normalization of ip is hxed by a physical 

reqnirement, it follows that dip is entirely determined by the shape of the membrane. 
28 

We conclnde that, to leading order, ineqnivalent ansatz are parametrized by the 
shape of the membrane and O on the membrane. 


A.5 Details concerning Myers-Perry black holes 

The metric for a rotating black hole in arbitrary dimensions is well known [23, 24]. 
The black hole metric is given in Eddington-Finkelstein coordinates by in 2N + 2 
dimensions 


N 


N 


ds^ = 2 idv — ^ttinldOi I dp + + a‘f){dpf + Pid6^) — dv^ 

\ i=l ) i=l 

N 

+ (^) {dv - ^ aipjdOiY + p'^da^, 


(A.ll) 


2 = 1 


while in 2N + 1 dimensions the metric is given by 


N 


N 


2=1 


ds^ = 2 ^ aipfdOij dp + + CLi){dp‘^ + PidO^) — dv^ 

N 

+ ip^i) - Pz 


(A.12) 


2 = 1 


In both (A.ll) and (A.12), the black hole is characterized by a mass parameter m 
and the N independent angnlar velocities Oj in the N angnlar directions 6i. Apart 
from 6i, the coordinates in both the metrics above consist of p, v the N ‘direction 
cosines’ pi yielding of total of 2N + 2 coordinates in all. Of conrse this is one 
coordinate too many for the odd dimensional metric, in which case the coordinates 

^®The fact that the normalization of dtp is physical is closely related to the following observation. 
dtp is proportional to dB. Under the rescaling B —^ aB, dB —^ adB on the membrane. However 
dtp is unaffected by this scaling on the membrane. 
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/ij are constrained to obey 


N 




(A. 13) 


In the case of the even dimensional black hole, the additional ‘coordinate’ a in (A.11) 
is given in terms of direction cosines by 

N 


a 




(A. 14) 


In both metrics the quantities 11 and F are dehned by 

N 

n = 


2=1 


F = 1 


N 

E 


2,,2 


(A.15) 


a^/i: 


^ + ^2 


(A.16) 


The rotating black holes can be recast in ‘flat space’ type coordinates as follows. 
In the case of the even dimensional black hole we move to the new coordinates 
Xi, Hi, z,t {i = 1... N) dehned by 

Xi = cos 9i — Qi sin 6 *j), 

Vi = pii{psin9i + aiCos9i), 
z = pa, 
t = V — p 

In the case of the odd dimensional black hole we use use the coordinates Xi, Ui and 
t dehned in a similar manner by 

Xi = pi{p cos 9i - at sin 9i), 

yi = Pi{psin 9i + ai cos 9i), (A. 17) 

t = V — p 

In terms of the new coordinates the even dimensional black hole A. 11 takes the form 

N 

ds^ = —dt^ + + dyl) + dz^ 

mp 


nF 




“ p^ + K / ^ p^ + a. 

4=1 ' / 4=1 

whereas the odd dimensional black hole A. 12 turns into 

N 

ds^ = —dt^ + ^(dx^ + dyl) 


- Vidxi) 


mp 

nF 


2=1 

dt+ll-f; 


N 


(A,19) 


“ p-* + a" 

4=1 ‘ 


dp-Y, 


2 I ^2 VidXi 


^ P^ + af 
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The quantities /i* and p in (A. 18) and (A. 19) are now functions of our new co¬ 
ordinates. In the even dimensional case these qantities are determined by by the 
equations 



(A.20) 


while in the odd dimensional case the relevant relations are 




(A.21) 


In this paper we are interested in the special case of rotating black holes with non¬ 
zero rotations turned on in only q two planes with q D. Given such a black hole 
(as mentioned in the main text) we choose p = 2q — 1. Let us now change notation 
so that Xi and pi are coordinates of the two planes with non-zero rotations, i.e. 
i = 1.. .q. We group the other Xi,yi coordinates (the coordinates of the remaining 
two planes in which the angular velocities all vanish) together as zm coordinates. In 
the even dimensional context (A. 18) these zm coordinates appear on equal footing 
with the coordinate 2 ; and we group that coordinate together with zm- With these 
dehnitions the index M runs over d + 1 values and the black hole metric, in both 
odd and even dimensions, takes the form 


<? 


ds^ = -dt^ + '^{dx‘1 + dy'^) + ^ dzl^ 



2 


(A.22) 


where /i^ = as above, and the equation that implicitly determines p is 



(A.23) 


where main text. 
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B Details of the perturbative procedure 

We move to a patch centred around Xq and use the following patch coordinates 


^ - 1), ^ - ij). 

^ = Ti(x„)(x>‘ - ij) 

In these coordinates the rescaled zeroth order metric G 
the following form. 

ds^ = 2soe-’^°dR dV - e-^’^°F{R)dV^ 


= S - So 

d 


(B.l) 


ab and the one form x fake 


X 



+ J2dyW + 0 

i=l 



^ + 2dR + 0 
So 





(B.2) 


where 

F{R) = 1^(1 - c^) - 

Ignoring corrections of order (B.2) is the black brane metric. 

In these coordinates the terms listed in the second line onwards of (3.8) takes 
the form 


^-R 


ds'^ 


Hi) 


Ki{R)dV^ + 2K2{R)dV dY + K3{R)dY^ 
2K4{R)dR dY + K^{R)dy^dy^ 

2Q\{R)dV dy^ + 2Q\{R)dY dy^ + %{R)dy^dy^ 


(B.3) 


As we have explained in the main text, at (^) the scalar, vector and the tensor 
data are effectively constants (equal to their values at the point Xq) as derivatives of 


^®(B.2) can be cast into the standard form (2.10) by the rescaling V = Y and shifting 

R = R + [2/io — log(l — (?)]■ The ‘a’ parameter of the standard form metric may be verified to be 
given by 

l-F 


In our computations below, however, we do not perform this scaling and shift; we work with the V 
rather than the V coordinate. 
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these expressions contribute to the metric only at O and higher. It follows that 


7 




A=0 

5 




71=1 

2 


T„(R) ^Y.Ta(R)T'M) 


7l=l 


Substituting (B.3) in Einstein equation and the scalar held equation as given in 
(2.2) we obtain ordinary differential equations for various metric functions. SO{p) 
symmetry ensures that the equations for the scalar functions Ka, vector functions 
Qa{R) and tensor functions Tij{R) are mutually decoupled. Of course the equations 
for all functions of a given symmetry (e.g. scalars) all mix. Nonetheless it turns out 
that the equations are rather simple and explicitly solvable 

Recall that our data included an antisymmetric tensor; as there is no antisym¬ 
metric tensor in the metric, it follows even without a calculation that there is no 
term in the correction metric proportional to this piece of data. 

B.l Scalar sector 

In this sector the relevant equations are Efiv, E^y, Ew, Eyy, Eyy and 

Etr = Here by E^^, we denote the second equation in (2.2). 

In order to solve these equations we found it useful to hrst study a linear com¬ 
bination of Erv and Eyv which turns out to give a decoupled equation for K 2 {R). 


El = soe’^°+^Evv + e^E{R)ERv 



A=0 


(B.4) 


where ' denotes the derivative w.r.t. R. The linear combination of equations pre¬ 
sented in the hrst line (B.4) evaluates (by explicit computation) to the second line of 
that equation. The hrst term in (B.4) is homogeneous in the unknown K^i^R) while 
the second term, & 2 {R), is the source term proportional to the scalar data S^. The 


the equations for the unknowns are all linear it is possible in principle- and useful in practice 
- to obtain solutions for pieces of one derivative background data turned on one at a time and then 
to add the results. 
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explicit expressions for these source terms are 




2 so 


R 


&\{R) 

&1{R) 



/ c®-6c^ + 8c2-1 

V 2 c 2 (c 2 - 1 ) 

^ 2 c^ H—^ — 6 — -R + 


{c^-2)R {c^ - 2) 

^ 2 (c2 - 1) 2 (c2 - 1) } 

(2c^ - 

J 


Gl{R) = ^ (1 - 

6 ^(R) = —, 6^(R) = e"°(-2 + c^) 

6 ®(R) = -2&1{R) = (1 - c^)e^° 


(B.4) is easily solved by integration and we obtain 


K^iR) 


2soe^F{R) 



(e ''&2{p)\ 

V npf ) 


(B.5) 


(B.6) 


The integration constant is hxed using the boundary condition that for large R, 
K 2 {R) should behave at most like a polynomial in R. 

Now the function R(R) has a zero at Rq = 2ho — log(l — c^). It follows that the 
integrand in (B.6) can be expanded around p = Rq as 


(e ^G>^{p)\ 

V np? ) 


Bt , 

(p-Ro)2 + p-Ro 


+ regular 


Upon integrating, it follows that 


R2(R) OC (R - Ro) 

A 


f BfSA 
\R — Rq 


+ B^Sa ln(R — Ro) + regular 


While the terms proportional to above are perfectly analytic, the terms propor¬ 
tional to B 2 are potentially singular (non analytic). It follows that K 2 {R) solutions 
is regular everywhere outside its event horizon (as we demand) if and only if 

J2b^Sa = o. 

A 


This equation may be rewritten as 

7 

Equation of motion oc E ^^(a«©^(R)|^=«J =0 (B.7) 

A=0 


upon explicitly evaluating the derivatives in (B.7) we hnd the scalar equation in 
(1.7). After imposing the equation of motion K 2 {R) is perfectly regular. 
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Next we solve E^y = 0. If we substitute the solution for K 2 {R) and the equation 
of motion, Ery becomes very simple: 

Ery = ^dn (e-^KsiR)) = 0 (B.8) 

/So 

The only normalizable solution to this equation 

KsiR) = 0 

Etr gives another decoupled source free equation for K^{R). 

Etr = (^e^E{R) dn (e-^MR)) ^ = 0 (B.9) 

If we impose regularity at R = Rq and normalizability at i? = cxd the only allowed 
solution is 

K,{R) = 0 

After we substitute the solutions for K 2 {R), KsiR) and K^i^R), the equation 
Err becomes the simple source free equation for K^i^R). 

Ery = ^dR {e-^K^{R)) = 0 (B.IO) 

/So 

From here it follows 

K^iR) = 0 

Finally we solve for Ki{R). We use the scalar held equation (i.e., the hrst 
equation in (2.2)) for this purpose. After substituting in the solution for all the 
other Ki{R) as well as the equation of motion, this scalar equation takes the form 


where 


7 

E^ = Y,Sa {dRKtiR) + &^{R)) = 0 

A=0 


(B.ll) 


<5\{R) 

<5\{R) 

&\{R) 

&\{R) 

&\{R) 

G\{R) 


-c® + llc^ - 23 c 2 + 11 + (c^ - 1)^ (3 c2 
(c2 - 2f (c2 - 1) 

0 


2so(l -2c^ + c2(i? + 3)) 

c2 (c"^ — 3c2 + 2) 

So (4c^-c"(/? + 7)-l) 

c 2 (c 2 - 2f 


So 

1 -c 2 ’ 
&1{R) = 


Gl{R) 

0 



1 


8 )i? 


(B.12) 
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(We have used the equation of motion to solve for Si in terms of the other scalar 
data. That is why = 0). 

Equation (B.ll) is trivially solved by integration and we hnd 

Ay(fl) = - r'dpef(p). 

Jo 

We have chosen the integration constant above to satisfy the arbitrary condition 
Kf{R = 0) = 0. A different choice of integration constant could be absorbed by a 
shift of "0 by a function of order J and so is physically irrelevant. 


B.2 Vector sector 


flere the relevant equations are Ejiyi, Eyyi and Eyyi- We hrst study the linear 
combination of Ejiyi and Eyyi which yields a decoupled equation for Q^^R). 

Ei = 2 (^e^°SQEyyi + soE{R)Ejiyi^ 

= J2SAUE{R)dn{e-^Q^{R)) +^t{R)] =0 

A=1 ^ ^ 


where 

= -so(2 - 

^1{R) = So(l - (B.14) 

= 0 , ml{R) = 0 , ^1{R) = 0 


Solving this equation with the boundary condition of ‘normalizability’ at inhnity we 
hnd 


pR ra 

Q2iR) = - / 
Jr 


dp 


vi(p) 

Fip) 


(B.15) 


Clearly Q^^R) above has a logarithmic singularity aX R = Rq] the coefficient of this 
singularity is proportional to the vector equation of motion 


5 

Equation of Motion = Vj^^^iRo) = 0 (B.16) 

A=1 


Upon substitution, (B.16) reduces to the second of (1.7). Note, however, that all 
the nonzero source functions TT^(i?) are proportional to each other as functions of R 
(they are all proportional to e~^). It follows that (B.16) implies that 

ys’opR) = 0 

A=1 


at all R and so it follows from (B.15) that 

Q2{R)=0 



Upon substituting the solution for Q 2 {R) and also the equation of motion, the equa¬ 
tion Ejiyi allows us to solve for Qi{R). It turns out that only Q{(-R) has a nonzero 
source and we hnd 

pho—R / \ p2ho—R 

= 0 (B,17) 

(B.17) is easily solved by two integrations. The solution is automatically regular 
everywhere. There are two integration constants. One of these is hxed by the 
requirement of ‘normalizability’ at i? = oo. The other integration constant can 
be absorbed into a shift in velocity held by a term of order We have arbitrarily 
hxed it by demanding that Qi(0) = 0. 

B.3 Tensor sector 

The equations are simplest in the symmetric tensor sector. The unique tensor equa¬ 
tion, Eyiyj, evaluates to 

EV^y^ = -J2[TAfy^e-^dR(e^E{R)dR {e-^TA{R)) ) = 0 (B.18) 

Quite remarkably no tensor source appears on the RHS of (B.18). The most gen¬ 
eral solution of this equation is easily determined, and turns out to be either non- 
normalizable or singular at i? = i?o- It follows that 

Ta^R) = 0 

and there are no corrections to the metric in the tensor sector. 
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